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New QCD sum rules for nucleons in nuclear matter are derived from a mixed cor- 
relation function of spin-1/2 and spin-3/2 interpolating fields. These sum rules allow a 
determination of the scalar self-energy of the nucleon independent of the poorly known 
four-quark condensates. An analysis of these new sum rules in concert with previous 
nucleon sum rules from spin-1/2 interpolators indicates consistency with the expecta- 
tions of relativistic mean-field phenomenology. We find M* = 0.68 ± 0.08 GeV and 
H v = 0.31 ± 0.06 GeV at nuclear matter saturation density. 
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m : 1 Introduction 

Understanding the observed properties of hadrons and nuclei from quantum chromodynam- 
ics (QCD) is a principal goal of nuclear theorists. The QCD sum-rule approach [|TJ is a 
particularly useful method for connecting the properties of QCD to observed nuclear phe- 
nomena. Recent progress in understanding the origin of the large and canceling isoscalar 
Lorentz-scalar and -vector self-energies for propagating nucleons in nuclear matter has been 
made via the analysis of QCD sum rules generalized to finite nucleon density M. These large 
self-energies are central to the success of relativistic nuclear phenomenology. 

However, the previous sum-rule predictions for the scalar self-energy are sensitive to the 
density dependence of certain "scalar-scalar" dimension-six four-quark condensates. This 
density dependence is unknown ||]. There are at least two ways to clarify the situation. One 
direction is to attempt to better determine the density dependence of the four-quark con- 
densates via modeling ||. An alternative approach, which is adopted here, is to derive new 
QCD sum rules for the scalar self-energy that do not depend on the four-quark condensates. 

The new sum rules are obtained from a mixed correlator of generalized spin-1/2 and spin- 
3/2 interpolating fields. The spin-1/2 states remain projected, and one generates additional 
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sum rules for the scalar self-energy that are independent of the problematic four-quark 
condensates. These new sum rules, along with previous sum rules from spin-1/2 interpolators 
that are also independent of four-quark condensates @, provide an interesting forum for 
examining the expectations of phenomenological approaches. In this work we present results 
from the first analysis of the new sum rules; more detailed analyses will be given in future 
publications j|. 



2 New QCD Sum Rules 



The finite-density QCD sum-rule approach focuses on a correlation function of interpolating 
fields carrying the quantum numbers of the hadron of interest 
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Here, the ground state of nuclear matter \ty ) is characterized by the rest-frame nucleon 
density and by the four- velocity . The nucleon interpolating fields are [^|, ^[ 



xi(z) = l^ ahc {[u T a {x)C-i b d b {x)] u c + l3 [u T a (x) C d b (x)] j 5 u c (x)} 



(2) 



abc 



{ ul(x)C p \d b (x) 



<y px iv u c {x) - 



v£(x) C a p \u b (x) a pX j v d c (x)} , (3) 



where T denotes a transpose in Dirac space, C is the charge conjugation matrix, and (3 is a 
parameter allowing for arbitrary mixing of the two independent spin-1/2 interpolators. For 
(3 — — 1, the correlator of \]i an d X v gi ves the sum rules discussed in Ref. 0. 

In the rest frame of the medium, the analytic properties of the invariant functions can 
be studied through Lehman representations in energy. The quasi-nucleon excitations are 
characterized by the discontinuities of the invariant functions across the real axis, which are 
used to identify the on-shell self-energies. A representation of the correlation function can 
be obtained by introducing a simple phenomenological Ansatz for these spectral densities. 

On the other hand, the correlation function can be evaluated at large space-like mo- 
menta using an operator product expansion (OPE). This expansion requires knowledge of 
QCD Lagrangian parameters and finite-density quark and gluon matrix elements. Finite- 
density QCD sum rules, which relate the nucleon self-energies in the nuclear medium to these 
QCD inputs, are obtained by equating the two different representations using appropriately 
weighted integrals 0. 

The correlator Uj^(q) contains nine distinct structures at finite density. Two of these 
sum rules are independent of the problematic "scalar-scalar" four-quark condensates and are 
dominated by the leading terms of the OPE. 

To analyze the sum rules, we follow the techniques introduced in Ref. || |(| for deter- 
mining the valid Borel region. We limit the continuum model contributions to 50% of the 
phenomenology, and maintain the contributions of the highest dimensional operators in the 
OPE to less than 10% of the sum of OPE terms. This defines a region in Borel parameter 
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space where a sum rule should be valid. Uncertainties are estimated via a new Monte-Carlo 
error analysis introduced in Ref. ||. 

As in previous works on finite-density sum rules, we use the linear density approximation 
for estimating the in-medium condensates. Central values and uncertainties of QCD vac- 
uum parameters are taken from Ref. 0. Finite-density parameters are in accord with the 
estimates of Ref. 0. Since the mixed condensate (g c qcr ■ Gq) PN is chirally odd, we assume 
the same density dependence as for the quark condensate. 

The one firm conclusion from previous in-medium studies is that the vector self-energy 
is positive and a few hundred MeV. Hence, we begin by fixing E„ = 0.25 GeV at saturation 
density and searching for a region in which both sides of the QCD sum rules are valid. 

There are two sum rules obtained from the correlator of \]i an d xl t na t are independent 
of the problematic four-quark condensates. One of these satisfies the validity criteria. 

+ 5 (5 + ^ 2 + ^ E q M 2 E (qUD oq ) PN L~ 4/9 
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+ g^2 q2 ^ iD o iD oq)p N + j^(9cq i cr-Sq) PN 

+ 5 + 2 ^ 5(32 E q (Jq) 2 pN L^. (4) 
Here Ai denotes the coupling of x\ t° ^ ne quasi-nucleon state. We have also defined 



M* N = M N + E S , E q = E v + y/q 2 + M* N 2 , E q = E v - ^/q 2 + M* N 2 , (5) 

where S s and £„ are the scalar and vector self-energies of the nucleon in nuclear matter, 
respectively. The anomalous dimensions of various operators have been taken into account 
through the factor L = ln(M 2 /A 2 ^ CD )/ln(/i 2 /A 2 1CD ) 0. We have also defined E = 1 - 
e -s /M 2 an( ^ jt; i = i _ e -s /M 2 ( SQ /j\/f2 _|_ w hich account for excited-state contributions [0. 

Of the two favorable new sum rules obtained from the overlap of spin-1/2 and spin-3/2 
interpolators, only the sum rule at the structure 7^^, satisfies the validity criteria. 

^W^ 1 '^ = -^^) PN M"E,L^-^{g c 7la.Gq) PN L-^ 



3 ' (qq} PN (^^) PN ^-L^. (6) 



96 ^™ 7i ,PN M 2 
The parameter (3 is selected to minimize continuum model contributions while maintain- 
ing reasonable higher- dimension operator contributions such that the pole may be resolved 
from the continuum contributions @. The optimal choice for the sum rule of @j is f3 = —0.7. 
The continuum model contribution of @ is independent of (3. Since the second term of (@) 
has little overlap with the ground state nucleon |7|], we set (3 = 0. 
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3 Implications 



In vacuum, higher-dimensional condensates are typically evaluated by adopting a factoriza- 
tion prescription. It is generally accepted that factorization fails for the four-quark operator. 
However, the dimension seven operator appearing in (§) appears to be reasonably well ap- 
proximated by the factorization prescription in vacuum |J. As a chirally odd operator its 
density dependence should be qualitatively similar to that of the quark condensate. Fortu- 
nately, this behavior arises naturally in the factorized form where the density dependence of 
the gluon condensate is estimated to be a 7% effect. Hence the large uncertainties associated 
with the density dependence of factorized operators are largely eliminated in these new sum 
rules. We will further examine the sensitivity of these assumptions in a later work M. 

The density dependence of the new sum rule of (|6|) is predominantly governed by the 
quark condensate and is common to all terms of the OPE. As such, the effects of increasing 
density will be to reduce the residue of the pole while the pole position remains largely 
unchanged. This result is the key feature absent in the former in-medium nucleon analysis 
0. The approximate invariance of M* + H v is manifest in (|6|). 

Fig. la displays the valid Borel regimes for the two sum rules (§) and @. The dis- 
tributions for M* and T, v are illustrated in Fig. lb. We find M* = 0.68 ± 0.08 GeV and 
T, v = 0.31 ±0.06 GeV at nuclear matter saturation density. These results are consistent with 
the expectations of relativistic phenomenology. 
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Figure 1: (a) Illustration of the valid Borel regimes for the sum rules of (|j) (large dash) and (|6|) 
(fine dash). Both continuum model contributions (limited to 50%) and highest dimension operator 
contributions (limited to 10%) are illustrated, (b) Histogram for £„ and M* obtained from a Monte 
Carlo sample of 1000 QCD parameters. 
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